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Abstract 

Stress analysis and preliminary design/optimization procedures are presented for gyroscope 
momentum wheel rings composed of metallic, metal matrix composite, and polymer matrix composite 
materials. The design of these components involves simultaneously minimizing both true part volume and 
mass, while maximizing angular momentum. The stress analysis results are combined with an anisotropic 
failure criterion to formulate a new sizing procedure that provides considerable insight into the design of 
gyroscope momentum wheel ring components. Results compare the performance of two optimized 
metallic designs, an optimized SiC/Ti composite design, and an optimized graphite/epoxy composite 
design. The graphite/epoxy design appears to be far superior to the competitors considered unless a much 
greater premium is placed on volume efficiency compared to mass efficiency. 

I. Introduction 

Gyroscopes use the conservation of angular momentum in order to measure or control orientation. 
Often referred to as control momentum gyroscopes (CMGs), these systems are used in aircraft, satellites, 
spacecraft, and ships (refs. 1 and 2) (see fig. 1) when they are cost and mass effective. A key component 
of a gyroscope is the momentum wheel, also sometimes referred to as the CMG flywheel. It is this 
spinning component that stores the angular momentum that can be used by the system for control or 
sensing purposes and thus, along with its hub, drives the design of the gyroscope system. While the non- 
terrestrial applications for gyroscopes demand an optimized light-weight and low-volume design, the 
momentum wheel design has not typically been optimized. For example, heritage designs have 
predominantly used a standard operating rotational speed of 6000 rpm on which legacy bearing designs 
have been based, despite the fact that higher speeds will provide significantly improved mass and volume 
efficiency (ref. 1). Designs are also often based on one material chosen a priori and a finite element stress 
analysis to determine the minimum margin under operating conditions (refs. 1 and 3). Clearly, this 
approach tends to be overly conservative, and a good deal of efficiency that could be captured through 
design and material selection is being left on the table. Further, with the development of efficient and 
reliable gyroscopes, additional system weight saving are possible by utilizing the gyroscope for energy 
storage as well as control (refs. 4 and 5). 

This paper presents an analytical stress analysis, applicable to both composite (anisotropic) and 
metallic (isotropic) gyroscope momentum wheel rings. The stress analysis is combined with an 
anisotropic failure criterion to enable the failure (rupture) prediction of the momentum wheel ring due to 
angular velocity and gimbal maneuver loading. A factor of safety is incorporated, and a sizing 
(optimization) procedure is presented, which was implemented in a computer code. While the method 
considers only a single rotating ring (rather than multiple momentum wheel rings, the hub, and the shaft), 
it is highly efficient and thus well suited for preliminary design and sizing studies in which hundreds or 
thousands of potential designs may be considered. Full system design, requiring a detailed finite element 
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(a) (b) 

Figure 1. — Example packaged CMGs from (a) the Space Shuttle Orbiter and (b) the International 
Space Station (mockup). 


stress analysis (cf. ref. 6), would be necessary to implement the designs considered herein into an actual 
gyroscope design. In the detailed design process, trades between the ring and hub would be considered to 
arrive at the overall optimum configuration. 

Results are presented for two metallic ring designs (stainless steel and AerMet 1 00) and two 
composite ring designs (SiC/Ti and graphite/epoxy) with unidirectional circumferential reinforcement. 
The results illustrate the competing nature of the need to minimize both mass and volume while 
maximizing angular momentum. A performance index, that combines these competing needs, is 
developed and applied to the design. Based on this performance index, the graphite/epoxy design emerges 
as the clear choice among the materials considered when an approximately equal importance is placed on 
both volume and mass efficiency. 


II. Stress Analysis 

An annular, circumferentially reinforced, gyroscope rotor momentum wheel ring that is considered is 
shown in figure 2. As shown, the ring may be subjected to internal and external pressure (possibly 
associated with press-fitting with additional rings), an angular velocity, co, and a linear temperature 
gradient defined by, 


A T = T a -T ref + {r-a)^-^- (1) 

b-a 

where a, b, T a , and T h , are the inner and outer radii and boundary temperatures (see fig. 2) and T ri , f is the 
reference temperature (refs. 7 and 8). We also consider loading due to an out of plane angle change at a 
rate of 0 (e.g., caused by a gimbal maneuver), as shown in figure 3. It is assumed that the out of plane 
stress caused by this gimbal maneuver is completely decoupled from the momentum wheel ring in plane 
response. Thi s assumption is valid provided the out of plane stress is much smaller than the in plane 
stresses, which, as will be shown, is true in the present application. The stress field that results from the 
problem posed above involves only normal stresses, with in plane components, a,, and c 0 , a function of 
radial position, r, and the decoupled out of plane component, a : , a function of distance from the rotor in 
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Figure 2. — An annular gyroscope rotor momentum 
wheel ring with circumferential reinforcement. 



Figure 3. — A gyroscope momentum wheel 
ring subjected to an out of plane angular 
change due to a gimbal maneuver. 

plane midline (e.g., a function of the coordinate y in figure 3). Under these conditions, the equilibrium 
equations reduce to, 


da r ! u , -g B 

dr r 


+ pco 2 r = 0 


which allows the circumferential stress to be expressed as, 


U / \ 2 2 

<*°=— (raj + pco r ~ 


( 2 ) 


( 3 ) 


Combining the equilibrium, compatibility, and constitutive equations, along with the assumption that the 
out of plane stress, o z , is small relative to the in plane stresses, the governing equation for the radial stress 
can be expressed as (see ref. 7 for details), 
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(4) 


d 2 a„ „ da 


+ 3 r ar + (' -z’H=e(r) 


dr~ dr 


where I = ^E L jE T , E L and E T are the longitudinal and transverse elastic moduli. 


Q{ r )~b a { r ^L ( L T b ^{2(x L a r ) + a\^E L (T b ^ c /)( a i a r ) + rpco (3 + v z )J ^ 

- b \_ E ,.{ T a - T r ef ){v- L -a r ) + r 2 par (3 + v z )]} 

a L and a T are the longitudinal and transverse coefficients of thermal expansion, 17 is the longitudinal 
Poisson ratio, and p is the density. The general solution of equation (4) is, 

T = Q'"' 1 + Cy h + Qi {/) - Qi [r) (6) 

and from equation (3), the circumferential stress distribution is, 

g 9 = Cj [m { + 1 )r"‘' + C 2 + l)r'" 2 + ( 777 , + 1 'jQ l (r) - ( 777 , + \)Q 2 (/) + pcoV 2 (7) 


where, 


m { 

Q\ ( r ) = — [ r^ m,+V] Q(r)dr 

m 2 

Q \ { r ) = — \r [m 2 + ' ] Q{r)dr 


( 8 ) 

(9) 


777 , = — 1 + E , 777 2 = — 1 — E , and C\ and C 2 are constants that can be determined from the boundary 
conditions. 

Limiting attention to the isothermal case with zero internal and external pressure, the normalized in- 
plane stress fields can be written as, 


g = 


pc o 2 b 2 


3 + v L 

S 2 -9 


l-X 


3+Et 


(X 2E -1) 


f-1 

yb. 


2-1 X 


X 2S_ X 3+S | f r 7 


-1-2 X 


(X 22 -l) 




+ 


f-1 

yb; 


f = 


3 + u, 




pc/Zi 2 I 2 -9 i^ 21 - 1 ) [b 


( [x^-X^Af rX^A 


(X 2S -1) 


+ 


(£ 2 +3u i ) ( 7- 
3+u, 


.X 2 


( 10 ) 


( 11 ) 


where X= alb and the normalized circumferential and radial stresses,/ and g, have been defined. In the 
case of isotropic materials, equations (10) and (1 1) reduce to, 


g = 


pwV 


3 + o 
8 


( fi-x 4 


X 


[(X 2 -1)J 


'(x 2 -x 4 V r Y 2> 



/ \ 

2 1 


r 

1 

+ 

— 

1 

( 


UJ 

! J 


( 12 ) 
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III. Failure Criterion 

Because the momentum wheel ring may be circumferentially reinforced, and thus transversely 
isotropic, a transversely isotropic failure criterion has been adopted that accounts for the multiaxial stress 
state. This multiaxial criterion is based on an equivalent transversely isotropic effective ff-likc) stress 
measure (multiaxial) that is comprised of physically meaningful invariants that represent stress states that 
are likely to influence strongly the various damage modes within composites (ref. 9). For example: i) the 
transverse shear stress (7)) accounts for matrix cracking, ii) the longitudinal shear stress { 1 2 ) dictates 
interfacial degradation, and iii) the maximum normal stress (/ 3 ) in the fiber direction dictates fiber 
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breakage. Note that these are not the standard invariants of the stress tensor. All three of the invariants are 
combined into a function representing the effective transversely isotropic A-Iikc invariant, 


where, 


F () = 


— ( 
V M 


(4^ ( ) 2 -l)/ 1 +- 






9 

4 


-^2 + . h 


I x — J 2 I + 2 3 

h=i-h 

/ 3 =/ 2 


S=- 


\ \ 


./, = -5,5, 


i=D v s jk s u 

Dj = d,dj 


S = a — a,, 5 .. 

ij ij ^ kk ij 


( 20 ) 


( 21 ) 


and di (i= 1, 2, 3) are the components of a unit vector denoting the local fiber direction. Note that within 
( ) one can insert either ultimate or yield measures of material failure (or limit load). ct ( ( and ct ( , are 

thus the longitudinal and transverse limit loads. In case of isotropy F reduces to classical^ failure 
criterion. 


IV. Optimization Procedure 

Noting that the normalization of the stress components (eqs. (10) to (13) and (18)) eliminates the eo 
dependence for / and g, but not h, the multiaxial failure criterion, equation (20), can be written in terms of 
the normalized stress components as, 

F ( , = ±1 } [ g 2 - 2gh(w) + /7(co) 2 ] + gA(co) -gf- Ma) + f (22) 

a ( )L 


where the eo dependence of h has been indicated. We now assume as before that the out of plane stress (h ) 
is small compared to the in plane stress components, resulting in, 


F () = 


pc.) 2 /; 2 

G ( )L 


Z, 2 () g 2 -gf + f 2 


(23) 


which can be solved for the maximum allowable angular velocity, 


co 


2 

max 


t ~< allowable _ 
F ( ) 


( )L 


P b 2 


J%)g 2 -gf + f 2 


where the allowable magnitude of the failure criterion function is given by, 


( 24 ) 
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F 


( ) 


0.999 

FS 


(25) 


and FS is the desired factor of safety to be employed in the momentum wheel ring design. A numerator of 
0.999 is employed in this expression, rather than 1.0 in order to accommodate a small, but finite value of 
the out of plane stress. Also in equation (24), the [ ] operator refers to the maximum value of the 

bracketed term occurring within the momentum wheel ring, which can potentially occur at any radial 
position. 

The momentum wheel ring’s angular momentum at the maximum allowable angular velocity is given 
by, 

H = I? ® max ( b ) = |[l ■ -X 4 y P t^ x (b) (26) 

where the functional dependence of gw on the outer diameter, b, has been denoted. Rearranging 
equation (26) results in, 


b 4 


2 H 


(27) 


Equation (27) can be solved by direct iteration for the momentum wheel ring’s outer radius, b, which, for 
a given desired angular momentum ( H ), rim thickness ratio ( X = alb), out of plane thickness (t), and 
material, will be the minimum outer radius value for the ring that can provide the desired angular 
momentum. This will also be the lightest and smallest ring that can provide the desired angular 
momentum. 

The momentum wheel ring sizing/optimization procedure can then be summarized as the following 
steps: 

1 ) For a given H, t,X, 0 , material, factor of safety, and choice of ultimate or yield limit stress, 
assume a value (i.e., initial guess) for the momentum wheel ring outer radius, bo. 

2) For iteration i, solve for the stress field throughout the momentum rim (this can be done using a 
sufficient number of radial solution points) and determine: 




(28) 


Note that the radial location corresponding to the maximum value of the bracketed term in 
equation (28) establishes the controlling failure point in the momentum wheel ring, and this 
expression is a function of current iteration’s value for the outer radius, b,. 

1 ) Determine the maximum allowable angular velocity associated with b 

1 

2 

(29) 


2) Determine a new value for the outer radius, 


Pfy 2 


t? allowable _ 

_ F () g (^ 

fel)g 2 -gf + f 2 
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\_ 

4 


( 30 ) 


2 H_ 

|rc[l - Jr 4 ]p?<o max (^-) 


3) If b has not converged, set i = i+ 1 and go to step 2. 

4) Perform final check of full failure criterion throughout the momentum wheel ring, 

F ( } } [g 2 - 2gh(<o) + h((o ) 2 ] + gh( co) -gf- Mu) + ^ (3D 

a ( )l ^*3 

Violation of this full failure criterion indicates that the out of plane stress may be significant and may 
call into question the validity of the stress solution. In such a case, the magnitude of the out of plane stress 
should be examined carefully. 

This procedure sizes the ring for a given angular momentum, which is typically the specification to 
which gyroscope momentum wheel rings are designed and manufactured. The material parameters, out of 
plane thickness ( t ), and maximum gimbal rate must also be known. The ratio of inner and outer radii 
(X= alb ) may or may not be known. As such, a computer program was written to size momentum wheel 
rings based on the above procedure where the value of A is automatically varied from 0.01 to 0.99. Plots 
of various quantities versus Aresult and can be used to find the truly optimized momentum wheel ring 
solution, even across various materials. 


V. Results and Discussion 

The materials considered for the gyroscope momentum wheel ring in this design study are AerMet 
100 (a high-alloy steel), custom 455 stainless steel, a graphite/epoxy composite, and a SiC/Ti metal 
matrix composite. The material properties for these materials are given in table 1. It is noted that 
transverse strengths for the composite materials represent static values that have not accounted in anyway 
for fatigue. Were the ring subjected to a significant number of loading cycles (e.g., spin up and spin down 
cycles), an appropriate knockdown of these transverse values would be required to account for matrix 
fatigue cracking. Alternatively, a multi-ring press-fit design, that would preload the ring(s) into 
compression, could be employed to alleviate the transverse fatigue cracking problem. 


TABLE 1.— MATERIAL PROPERTIES. FOR COMPOSITE MATERIALS, DOUBLE VALUES REPRESENT 
LONGITUDINAL AND TRANSVERSE VALUES, RESPECTIVELY 


Material 

Volume 

fraction, 

(%) 

P 

(lbm/in. 3 ) 

Ultimate strength, 
(ksi) 

Yield strength, 
(ksi) 

Elastic modulus, 
(Msi) 

v L 

Stainless steel 

- 

0.284 

185 

175 

29.0 

0.26 

AerMet 100 

- 

0.285 

294.8 

247.4 

28.7 

0.30 

Gr/Ep 

60 

0.057 

302 

10 

302 

10 

23.1 

1.3 

0.282 

SiC/Ti 

35 

0.147 

201.6 

60.9 

201.6 

39.6 

26.1 

17.5 

0.272 


Throughout this design study, a factor of safety of 2 and a gimbal rate ( 9 ) of 1 rad/sec were assumed 
based on internal communications. Note that a separate study (not shown) indicated that in order for the 
out of plane stress to reach the magnitude of the material failure stresses, a gimbal rate on the order of 
1700 rad/sec would be required. This is clearly unrealistic, and in all cases with the 1 rad/sec gimbal rate, 
the out of plane stress was small, on the order of 0. 1 percent (at most) of the in plane stresses. Also, all 
results presented have employed the ultimate strength in the design. 
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We first consider the AerMet 100 material to illustrate how the results can be used to optimize the 
gyroscope momentum wheel ring design. Figures 4 and 5 show that, as one would expect, higher 
momentum wheel rings must be larger, and as the rim thickness decreases (higher X values), the outer 
radius of the ring must be larger to achieve a given momentum. 

A better assessment of momentum wheel ring’s size is its package volume, defined as, 

Package Volume = PkV - nb 2 t (32) 



Figure 4. — Optimum outer radius as a function of inner to outer radius ratio (X) for 


AerMet 100, t- 1 in., and various design-to angular momentum values. 



Figure 5. — Optimum outer radius as a function of inner to outer radius ratio ( X ) 
for AerMet 1 00, H = 1 700 ft-lb-sec, and various out of plane thicknesses. 


NASA/TM- 


-2007-214967 


9 


This is not the material volume of the momentum wheel ring, but rather the real volume occupied by the 
ring, where the open internal region contributes to the volume. This is important because it is this volume 
that must be designed for in housing the gyroscope. Figures 6 and 7 show the package volume of an AerMet 
100 momentum wheel ring. The fact that these curves are flat up to approximately X= 0.60 indicates that no 
further volume can be saved by reducing the internal to external radius ratio lower than this value. By the 
same token, raising this ratio beyond this value brings with it a package volume penalty. 



Figure 6. — Optimum package volume as a function of inner to outer radius ratio 
(X) for AerMet 100, f = 1 in., and various design-to angular momentum values. 



Figure 7. — Optimum package volume as a function of inner to outer radius ratio 
(X) for AerMet 1 00, H = 1 700 ft-lb-sec, and various out of plane thicknesses. 
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Figures 8 and 9 show the mass of an AerMet 100 momentum wheel ring. The trend is opposite that of 
the package volume curves, with higher X values resulting in lower masses. This is because higher X values 
are associated with lower rim thicknesses and higher outer radii. Large and thin rings provide the best 
angular momentum per unit mass (as the mass is located further from the center of the ring), but clearly, 
from figures 6 and 7, a trade off exists between mass and package volume. Both of these factors are 
important for design as higher mass of the ring is detrimental to the vehicle as is higher package volume due 
to additional mass required to house a larger gyroscope ring and accompanying system hardware. 



Figure 8. — Optimum mass as a function of inner to outer radius ratio ( X) for AerMet 
100, t = 1 in., and various design-to angular momentum values. Note that the 
value of b varies (as indicated) along each curve. 



x (= a/b) 

Figure 9. — Optimum mass as a function of inner to outer radius ratio ( X) for AerMet 
100, H = 1700 ft-lb-sec, and various out of plane thicknesses. Note that the value 
of b varies (as indicated) along each curve. 
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In order to quantify and trade the effects of mass and package volume, we consider the angular 
momentum provided by a ring per unit mass ( H/m ) and per unit package volume ( H/PkV) as quantities 
that should simultaneously be maximized. Plots of these two quantities vs. each other are given in 
figures 10 and 1 1 for AerMet 100. Because these plots are on a log-log scale, a line of slope 1 
corresponds to designs of equal package volume/mass ratios. Such lines also represent equal values of X. 


1000 


E 

n 


£ 100 



10 


H = 95 ft-lb-s 
H = 325 ft-lb-s 
H = 700 ft-lb-s 
H = 1700 ft-lb-s 


10 

H/PkV (in-lbf-s/in A 3) 


100 


Figure 10. — Angular momentum per unit mass plotted versus angular momentum per 
unit package volume for AerMet 100, t = 1 in., and various design-to angular 
momentum values. Data points in the upper right hand corner represent the best 
designs. 



Figure 1 1 . — Angular momentum per unit mass plotted versus angular momentum 
per unit package volume for AerMet 100, H = 1700 ft-lb-sec, and various out of 
plane thicknesses. Data points in the upper right hand corner represent the best 
designs. 
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Lines of slope -1 correspond to equal values of ( H/PkV) x ( H/m ). Because it is desirable to maximize 
both of these terms, this concept can be generalized to develop a performance index ( PI) for the 
momentum wheel ring, 


PI = 


f H ^ 

a 


{pkv) 


l m) 


(33) 


which, on a log-log plot like those provided in figures 10 and 11, yields lines of slope -a/p corresponding 
to lines of constant PI. The higher and further to the right that the line falls, the better the performance 
index and the better the design. Therefore, the relative weighting of the importance of package volume 
and mass minimization with respect to angular momentum can be quantified by the ratio of a to p. For 
instance, if package volume and mass are given equal importance to the momentum wheel ring design, a 
slope of — 1 results. As shown in figure 10, the shapes of the curves for the different angular momentum 
values indicate that an optimum design exists that maximizes this performance index. 

Figures 12 and 13 now plot this equally weighted performance index, equation (33) with a = P = 1, 
versus the inner to outer radius ratio, X. Clearly, in all cases an optimum design can be readily identified 
as the maximum in the curve. It occurs atX= 0.74. This is essentially the point at which the benefits of 
lower mass associated with higher X values (see figs. 8 and 9) overcome the benefits of lower package 
volume associated with lower X values (see figs. 6 and 7). 



x (= a/b) 


Figure 12. — Equally weighted performance index (see eq. (33)) as a function 
of inner to outer radius ratio (X) for AerMet 100, t = 1 in., and various values 
of design-to angular momentum values. 
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Figure 13. — Equally weighted performance index (see eq. (33)) as a function of 
inner to outer radius ratio ( X) for AerMet 100, H= 1700 ft-lb-sec, and various 
out of plane thicknesses. 


It is now possible to compare the optimum momentum wheel ring designs for the four materials 
whose properties are given in table 1 (Gr/Ep PMC, SiC/Ti MMC, AerMet 100, and Stainless Steel). 

Figure 14 shows a plot of the equally weighted performance index versus inner to outer radius ratio ( X) in 
the case of an H= 1700 ft-lbf-sec momentum wheel ring with an out of plane thic kn ess of 1 in. Recall 
that a factor of safety of 2 has been employed, and the momentum wheel rings have been sized to 
ultimate. This figure shows that the stainless steel design is inferior over the entire range of X values. At 
lower X values (below 0.4), the remaining three materials are comparable in terms of performance index. 
As the value of X increases, the PMC material’s performance index increases much more rapidly than 
those of the other materials and reaches a maximum value that is approximately 75 percent higher than 
that of the MMC and AerMet 100. This is due to the extreme low density of the PMC. Also, the PMC 
reaches its optimum design at X= 0.81, while the other three materials reach their optimum designs atX= 
0.74. The optimum PMC design corresponds to an outer radius of 16.7 in. and has a maximum tip speed 
of 844 m/s. 

Figure 1 5 provides a log-log plot of H per unit package volume vs. H per unit mass for all four 
materials and all four angular velocities considered previously and an out of plane thic kn ess of 1 in. 
Utilizing a line of slope — 1, which corresponds to the equally weighted performance index, it is clear from 
this figure (as was shown in figure 14) that the PMC material provides the best design (as measured by 
this performance index). Also plotted is a line of slope -2.7 that intersects an extreme point of both the 
PMC and AerMet 100 curves. This indicates that in order for the AerMet 100 material design to be 
competitive with the PMC material design, the importance of package volume must be 2.7 times the 
importance of mass to the design of the gyroscope momentum wheel ring (see eq. (33)). 


NASA/TM— 2007-2 14967 


14 


30000 



x (= a/b) 

Figure 14. — Comparison of the equally weighted performance index (see eq. (33)) 
as a function of inner to outer radius ratio ( X) for the Gr/Ep PMC, SiC/Ti MMC, 
AerMet 100, and Stainless Steel, t = 1 in., and H = 1700 ft-lbf-sec. 



Figure 15. — Comparison of angular momentum per unit mass plotted versus 
angular momentum per unit package volume for the Gr/Ep PMC, SiC/Ti MMC, 
AerMet 1 00, and Stainless Steel, f = 1 in, and various values of H. 
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VI. Conclusion 


A stress analysis and preliminary sizing/optimization procedure for monolithic and composite 
gyroscope momentum wheel rings has been developed. This enables improved design efficiency for these 
critical components that are found in many aerospace and terrestrial applications, and that have typically 
not been fully optimized in the past. The developed methods and results indicate the presence of a 
competition between the need to minimize mass and minimize volume, while maximizing angular 
momentum in both cases. A performance index was developed that combines and quantifies these 
competing design drivers, given the relative weighting of their importance. Results for two metallic 
(stainless steel and AerMet 100) and two composite (SiC/Ti and graphite/epoxy) designs indicated that, 
based on equal weighting of the volume and mass design drivers, the graphite/epoxy momentum wheel 
ring represents the most efficient design. This is mainly due to the high specific strength of the material 
compared to the other materials considered. If a much greater importance is placed on volume efficiency 
compared to mass efficiency, the AerMet 100 design becomes competitive, but only if the relative 
weighting factor for volume is greater than 2.7 times that for mass. 
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